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I. INTRODUCTION 



In the first part of this work (hereafter referred as M) we have established an exact field theoretical representation 



Q of the statistical mechanics of the restricted primitive model of electrolytes (RPM), i.e. an equimolar mixture of 
O , positively and negatively charged hard spheres of the same diameter a. Although the generalization to assymmetric 
models is certainly possible we have limited ourselves to the case where Ahe cations and the anions bear opposite 
I ' charges ±e. By applying the Kac-Hubbard-Stratonovich-Siegert-EdwardsoLl transform to the Coulomb interaction 
^ , we have obtained an exact new expression of the action S[4>\ rVffy similar to the well-known Sine-Gordon action to 
' which it reduces in the limit of vahishing hard-coi:e.diametersJ3'Ll 

] By contrast with recent works in that directionEJ^M our formalism incorporates correctly the hard-core effects and 
our Sine-Gordon action is well defined as a consequence of a smearing of the charge of the ions which regularizes the 
^ Coulomb interactions at short distances. In the companion paper we have applied our formalism to a derivation of 
1 the equation of statej-afjibe RPM in the low activity (z) regime. The old-results obtained in the frame of the Mayer 
' diagrammatic theorylljilJ and recently rediscovered by Netz and OrlandEl are confirmed and our derivation precises 
r the domain of validity of the (approximate) field theoretical formulation of Netz and Orland. 

In the present paper we show that the formalism developped in I also enables a derivation of the equation of state 
^ ^ of the RPM in the high-tempcratupe|Keffime and we recover the results obtained by Stell and LebowitzE3 in the frame 
I ' of the 7-ordering theory of liquids .llZlila The latter is a sophisticated diagrammatic perturbation theory in which the 
' O ' "small" parameter is 7 = /3e^ (/3 inverse temperature), otherwise the density is arbitrary. Note that a derivation of 
^ ■ the equation of state of the RPM at high density is clearly out of the scope of the theory of Netz and Orland in which 
hard-core effects are incorporated only at the level of the second virial coefficient. 

Our paper is organized as follows. In Sec. ^ we reorganize the Sine-Gordon action of the RPM as an expansioti. 
in powers of 7. A 7-ordering of the grand-potential lorpm is then obtained by making use of the cumulant thcorcm.tHI 
Note that we thus make use of the same technics as in I to obtain the ^-ordering of ojrpm which confers an aesthetic 



5_] • unity to our approach. In Sec. [II we derive an expression of wrpm valid at order 7^/^. Moreover we show that our 
expression for urpm is independent of the precise nature of the smearing (i.e. independent of the smearing diameter 
a, < a < (t) at this order. This result completes the independence of the z-expansion of lorpm upon a which was 
proved in the companion paper. Admitting this independence at higher orders in 7 greatly facilitates the algebra of 
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Sec. [V where lorpm is obtained at order 7^/^ by performing the calculations with a ^ 0. A Legendre transform 
then gives-an expression of the free energy of the RPM at order 7''/^ which coincides with the old result of Stell and 
Lebowitz.ll3 We conclude in Sec. ^ 



II. 7-ORDERING OF THE GRAND-POTENTIAL 



Our starting point will be the expression of the grand potential per unit volume ujrpm{P,i^) = 
— In [SflPM {P,t^,L^)] /L^ { (ii I/. L"^) grand-partition function) of the RPM, which was derived in I in the 

frame of a regularized Sine-GordonDfl#i'&i theory. Here i> = (3fj, (/? inverse temperature), where fj, = fj,± denotes the 
chemical potentials of the anions and the cations which can be chosen equal for simplicityu Moreover, we implicitly 
assume that the ions are enclosed in a cubic box of side L, with periodic boundary conditions. Recall that (c/. Eq 
(3.1) of ref. I). 



0.1 = - 



(27r)3 
ln(exp {~H 



In l + 7poW^r(9) 



(2.1a) 
(2.1b) 
(2.1c) 



In eq. (2.1a) ajHsii^a) is the gra nd po tential (per unit volume) of a fluid of neutral hard spheres at the chemical 
potential j/q = + ln(2). In Eq. ( |2.1b ) we have denoted by po = PHsi^o) = ^d^HS I ^vq the number density of this 
reference system. In Eq. (2.1b) Wriq) is the Fourier transform of the electrostatic interaction Wr{r) between two 
balls of diameter a distant of r, each of unit charge and with an uniform surfacic charge density I/tto^. The smearing 
of the charge over the surface of a sphere of diameter a smaller than the diameter cr of the ions was introduced in I 
in order to regularize the Sine-Gordon transform. Since 



(2.2) 



the expression ( 2.1b ) of the contribution wi is convergent for a 7^ as a consequence of the smearing. It must be 
stressed that the expression (2.1) of ojrpm is exact and does not depend on the value of the smearing diameter a (as 
long as < a < <t) 



The last contribution (2.1c) to the grand potential involves the effective hamiltonian 

'H[4>] = u[4>]-u„[<p] , 

C/[<^] ^ - ^ ^ / d^ri • • • dV-„ ft'") (1, ■■■,n)Y[ [exp(7/a) cos ( V7'/'(*)) - 1] , 



(2.3a) 
(2.3b) 

(2.3c) 



n=l 



i=l 



where 7 — (3e^ (±e charges of the ions) and where the h^\l, ■ ■ ■ ,n) are the connected n-body correlation functions of 
the reference hard sphere fluid at the density pq. These functions, which are closely related to the Ursell dis tribu tion 
functionsjl3 are solution of a hierarchy of integral equations which is briefly rewied in Appendix In Eq. (2.1c) the 
brackets (. . .)x^ denote a Gaussian average over the fluctuations of a real scalar field 4> with covariance X^-; the latter 
is a screened potential related to the Coulombic potential Wr by the relation 



Xriq) 



Wriq) 



l+JPoWriq) 



(2.4) 



In the companion paper I, we have given a low-fugacity expansion of lorpm, here we rather consider the behavior 
of the grand potential at high temperatures and arbitrary chemical potential v. In that aim we first^eorganize the 
Sine-Gordon Hamilton ian "H [(/>] as an expansion in the (small) parameter 7 of Stell and Lebowitalj. It is easily 
established from Eqs. (2.3) that 



(2.5) 



N=2 



2 



where 



N 



Wn 


n=l 


An) 


n\ J 




(-1)'* a"-P 
(2s)! {p-s)\ 



l2s 



(2.6a) 
(2.6b) 
(2.6c) 



The summations in the r.h.s of Eq. (2.6b) are restricted to all the integers pi> 1 {i — 1, . . . ,n) such that X]"=i Pi — 
which was denoted by the symbol ■ Later, we shall need explicit expressions of and W3; they are given in 
Appendix ^. 



In order to obtain a 7-expansion of Ao; we first apply the cumulant thcorcinc2l to Eq. (2.1c) which leads to 



ln(exp {~n m)x. = ^ - {n m)x..^ + 
a 2 



(2.7) 



and then re-order the terms in the r.h.s with the help of Eq. (2_^) which yields 



L^Aw = - ln(exp (-7^ [0]))^^ 



u 



Note that the subscript "c" in Eqs. ( |2/7| ) 

^n}x 



and {W^)x^,^ = {W^)x^ - {Wn)^,^. Gathering Eqs. (U 



|2.8j) denotes a cumulant average, for instance we have {Wn)x^_. 
UM) we can write finally 



UJRPM iP, v) = ojhs{vq) + 5^1 + Aw , 
Po7 
a 



Sloi = uji 



Y(W^2>x.,.+0(7') 



(2.8) 

{Wx)x. 

(2.9a) 
(2.9b) 

(2.9c) 



Eq. ( 2.9c ) needs some comments. The first order cumulants {WN)x^a which enter the equation originate from the 
7-ordering of (? Y[r/)]) Y, which is a regular function of 7 in the limit 7 — s- 0. Indeed {Ti.[(j)\)x^ ^ is obtained by applying 
Wick's theoremc^OEJ to Eq. (2^) and, since in the limit 7 ^ we have Xr Wr the resulting expression of 
{'H[4>\)x^a will involve integrals of products of long-range Coulombic terms by the n-body correlations /iq"'' of the 
hard sphere reference fluid which are assumed to be short-ranged and make the integrals convergent. Therefore the 
first order cumulants {Wn)x^ ^ regular in the limit 7 — > but it could be that some higher order cumulants are 
s ingu lar. For instance, it will be shown in Sec. |^ that c diverges as 7^^/^ in the limit 7^0. In writing Eq. 



(2.9c) we have assumed that, at this order, higher order cumulants do not contribute to AZU, which is, admittedly, 

difficult to prove rigorously. 

We stress again that, although the variable "a" appears explicitly in Eqs. (2.9a) as well as in the expressions ( |2.6| ) 
of the Wn, the grand-potential ujrpm does not depend upon "a" , Va e]0, cr]. Moreover, in a perturbation scheme such 
as the 7-ordering this independence should certainly be satisfied for each term of the expansion. We have checked this 
point up to order 7^/^. This result, which is discussed in details in the next section, completes the results of paper I 
where we have checked that each term of the low activity expansion of uirpm was actually independent of a. 



III. 7-EXPANSION OF lurpm{i3,u) UP TO ORDER 7 



i/2 



Since, in Eq. ( ^.9a ), ojhs{vo) does not dep end upon the temperature, we are left with the problem of finding a 
7-expansion of the terms 5uji and AZU of Eq. ( 2.9a ) up to order 0(7^/^). We have shown in I that, for a given a, uji 
was a function of the sole parameter (po7)^^^; therefore the expansions of oji in powers of (po7)^^^ given in I are valid 
either at low activities (or equivalently at low pg) or at high temperatures (i.e. at low 7), which is the case that we 
consider here. Recall that (cf. Eq. (3.6) of I) 
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(3.1) 



Note that SuJi is regular in the-Jimit a — > and that hnia^o Slui = —2 (npoj)'^^'^ /Stt which is reminiscent of the 

weh-known Debye-Hiickel resultESEJ. 

On the other hand, up to order 0(7^/^) we have, from Eq. ( 2.9cD 



L^AcJ = 7^(W^2)x.,. +0(7') 



(3.2) 



The expression of W2 is given in Appendix^. Taking the Gaussian average of Eq. (Bla) on finds 



2 



Po 
4! 



(3.3) 



where Pq is the derivative of the number density of t he hard sphere reference fluid po with respect to the chemical 
potential I'o. The Gaussian averages which enter Eq. ( |3.3| ) are computed with the help of Wick's theorem (details are 
given in Appendix]^). Inserting Eqs. (CI), (C2), (C5) in Eq. (jsj) yields 



where /iq^''(?'i2) = /iq'^''(12) is the usual pair correlation function in the limit of an infinite homogeneous and isotropic 

(2) I K 

system. Note that, in the core, i.e. for < r < a, {r) = —1. In order to derive Eq. (3.4) we have also made use 
of the compressibility sum rule for the reference hard sphere fluidcfl, i.e. 



.(2) 



(3.4) 



(3.5) 



- (2) (2) 
where (q) is the 31? Fourier transform of /iq (r). We note that the screened potential Xr (r) is a function which 



depends on the smearing diameter a and therefore, a priori, {W2)x^ as given by Eq. (3.4) depends upon a in a 
complicated manner. It is this dependance that we want to study in detail now. 

We have shown in the companion paper that, for a given smearing diameter a, the screened potential X-^ir) was a 
function of the sole variable 



Ko = 2(^po7)'^' 



(3.6) 



Note that kq is not, strictly speaking, the inverse Debye length of the system since it involves the density po of the 
reference system which is not equal, a priori, to the density p of the RPM. We have shown that, for r = (c/Eq. 
(3.8) of I) 



XAO) = -- 2(7rpo7)^/' + ^(^P07)« - ^i^Polf^'a' + Oij'a') 
a 15 3 



from which it follows that the first term of the r.h.s. of Eq. (3.4) is regular as a — > and reads as 

, 2 / 1 ^r(O) X^{Q)\ 1 ' 3 7/2 ^ 

For r 7^ 0, Xr{r) is a piecewise function defined as : 



X>(r) = q^{Koa/2)^^^^^-^ + 0{ja) for r > a , 
r 

X<ir) ^-(l-^)-Ko + 0(7a) for r < a , 



(3.7) 
(3.8) 

(3.9a) 
(3.9b) 



where q(x) = sinh(a;/2)/(x/2). Note by passing that Eqs. (3.9) imply that 



lim X-rir^ = y(r) = exp(— Kor)/r Vr. 
a— >o 



(3.10) 
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Inserting the expressions (3.9) of Xr{r) in the integral of the r.h.s. of Eq. ( |3.4[ ) and noting that h'^^ 
< r < a < a yields 



-1 for 



d^fhlf>{r)X>^r) 



(3.11) 



where the terms proportional to a and in the r.h.s. of Eq. ( [3.1l[ ) originate from the integration of X^'^{r) — X^'^{r) 
in the core (where h'if' (r) — —1). Moreover, note that, in Eq. ( |3.11 ) we can safely expand X^{r) in powers of kq since 
h'j^^ is a short range fmiction of r. Performing this expansion and gathering Eqs. ( |3.l| ), ( ^.^ ) and (3.11) we obtain 
finally, at order ^^^^ 

,(2) 



/a \ i \ ^ ( ^3/2 (Po7) 

Att 4 



(3.12) 



Due to a complex mechanism of compensations, the expression ( 3.12 ) of the grand potential at order 7^/^ is, as 
anticipated, indeed independent of the value of the smearing diameter a. 



IV. 7-EXPANSION OF THE GRAND POTENTIAL AND THE FREE ENERGY UP TO ORDER 7^^^ 
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In this section we obtain the expansion of u}rpm{P,v) up to order 7^/^. In order to make the derivations more 
tractable, we shall assume that the terms of order greater than 7^/^ in the 7-expansion of iORPM{f3,i^) are still 
independent of a. Moreover we shall consider the limit a ^ of Eq. (2.9a) which allows us to replace Xr{r) by its 
limiting expression y{r) = exp(— Ko?')/r and greatly facilitates the derivations. Therefore one has 

(4.1a) 
(4.1b) 



i^Fi.PAi iP, I') = lim (ujEsii^a) + SllIi + Acj) 



■0(7') 



The limit of Scui for a — > is easily obtained from Eq. (pT^) and reads as 



lim Sloi 



^=-i^i-poi) 



3/2 



(4.2) 



The expressio n of the term (W2)xt ^ the r.h.s of Eq. ( 4.1b ) has already been derived in Sec. [II. We take the limit 
a of Eq. (3.4) which gives 



lim7'(I^2)x_i"^ 



a— >0 



Wo7^_7^ /" ■3^,(2). x exp(-2Kor) 
2 4 J W ^2 

o2„,2 



+ 



87r3/2 



(po7)'^' / d^fhl^\r)r + Oi^'), 



(4.3) 



where we have expanded the Yukawa potential in powers of kq: taking advantage of the assumed short-range behaviour 

of hl^\r) and also made use of the compressibility sum rule ( |3.5| ). 

The expression of W3 is given in Appe ndix Taking the Gaussian average of Eq. (Bib) and then making use of 
Wick's theorem and of the relation ( AS! ) of Appendix W yields 



1 

"4^ 



pI 



J dVadVa 4''^(1,2,3)X,(1,2)X,(1,3)X,(2,3) 



lim w^ia) 



Po ^0 
48 



(4.4) 
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We take now the limit a — s- of Eq. (4.4) and expand the resulting expansion in powers of kq which gives 



a^o • 6 J ri2ri3r23 4 J dvo\ 



2 



(3) 

(P07)'/' / d?r,d^r, ^^^LiilM + o{^^) . (4.5) 



,(3) 
)_ 

Note that both cumulants {W2)xr ^nd (W^a)^^ are regular in the limit 7 — > which confirms the analysis made at 
the end of Sec. ||. 

The analysis of the limit a ^ of the last contribution —"/^{W2)x^ aK'^^'^) of the r.h.s of Eq. ( 4.1bD is more delicate 
and postponed to Appendix One finds that this term can be written as 

(POl)"^ f ^3^ ,3^ ,3^, (2), 



hm 7U-^(W^2^)x_ = ^7 V/2 1™ / rf'^2dV3d^f4 hl'\\,2)h\,'\iA) x 

O— >0 ' 47r ' Kq— ►0 J 

Ko exp {-Kq (ri2 + r23 + r34 + r4i)) 



''I2''23''34f41 



+ 0(7*). (4.6) 



By contrast with (W2)x-r and {W3)xr the cumulant {W2)xt ^ is singular as 7"^/^ in the limit 7^0. It is of course 
very difficult to give a precise analysis of the singularities of higher order cumulants. 

Gathering the intermediate results we get our final expression for the grand potential of the primitive model up to 
order 0(j'^''^). 

LORPAiiP, v) = uoHsivo) + Auj{P, v) + 0(7^) , (4.7a) 
9 

Au;(/3,j.) =^c.« , (4.7b) 

1=1 
2 



UJ 



UJ 



/dV^, (4.7d) 

(2)( 

r 



(3) = ^1/2 J KM , (4.7e) 

a;W = npo (^2po - ^p^^ 7' , (4.7f) 

J^^^i^[d^f2d^,,lil^, (4.7g) 
6 J ri2ri3r23 



,(6) = 

3 



= -(^7^/^ / d^r^(plh^^\r)) 1 , (4.7i) 



ri2r23 



= J-q^lunJ d'f2d'rsd'f, ^^^(l, 2)4^)(3, 4) x 



,(9) _ 

Kq exp (-Ko (ri2 + r23 + r34 + r4i)) 



''12?'23»'34?'41 



(4.7k) 



where we have made a change of notation and redefined Aui. Tinf free energy fRpuiPi p) is defined as the Legendre 



at 



transform of ujrpm{(3, v) with respect to the chemical potentials, i.e. 

fRPAliP, P) = OjRPMiP, v) + pif3, v)v , (4.8) 



where /5(/3, v) is the density of the RPM which can be obtained by differentiating —lorpm as given by Eq. (4.7a) with 
respect to vq. Therefore we have 
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(4.9) 



It should become clear in the sequel that an expansion of Ap to order 7^/-^ is sufBcient for our purpose. Let us denote 
by 

I/* = z/Q + Ai/ , (4.10) 
the chemical potential of a fluid of hard spheres at the density p. Eq. (^^) suggests that we search At/ under the form 

v3/2 ^ror^E)/^^ (4.11) 
In order to find the expressions of V'i/2 and V2 we note that, by definition of l^u one has 

P = PHsi^a + Ai/) 



v5/2^ 



Comparison of Eqs. ( [l.9| ) and ( 4.12 ) yields, by identification 



1/2 1/2 



1'2 = 



4Po 



(4.12) 

(4.13a) 
(4.13b) 



We have now in hand all the ingredients to compute the free energy. Making use of Eqs. (4.9), (4.10) and (4.11 
finds, at order 7'"/^ 

/fiPM(/3, p) = (^RPAiiP, v) + pv 

= ^Hs{v* — Aj/) + pv + Aw 
1 , 



i^Hs{^*) + pAi/ - -Po(Aiy)2 +pv + /^uj + 0(7-*) 
/h5(p) + /^'^ Ac. + 0(7*) , 



(4.14) 



where fnsip) = ^^Hsiv*) + p{v* — log 2) is the free energy of the reference hard sphere fluid at the same density p 
than that of the RPM (note that the free energy of the reference fluid and that of a fluid of identical hard spheres 
differ by a mixing entropy term) and 



/('^ = -^^^'2^^3/27'/' 

In order to obtain a more transparent expression of Jrpm we first remark that 



(2) 



1 



pg4^)(r)+.V2^1/2^3/2^(^2,(2)(^) 



(4.15a) 
(4.15b) 

(4.16) 



and that the term in brackets in the integral can be interpretated as a Taylor expansion of the function p'^h^^^\r). 
Since, from Eq. (4.11) Ai^ = 7r-^/2pJ/2^3/2 _|_ 0(^2-j^ rewrite (4.16) as 



^(2) + Jl^J = _ 



2 2 

(2) _ 7 P 



Oil') 



(4.17) 



where we have denoted by h*^^\r) the pair correlation function of the hard spheres at the density p (i.e. that of the 
RPM). 

Our second remark is that the combination of several terms occuring in Eq. (4.14) give the usual Debye-Hiickel free 
energy — K^/127r. Here n = 2(7rp7)^/2 is the true Debye length of the RPM which differs from kq by terms of order 
7'^/2_ Indeed one has 



7 



12^ 



^1/2^3/2^3/2 



3 

2_i 



3 po 



(4.18) 



By combining Eq. (4.1J) and the expressions of wi, uj^, uj^ one finds 



^(l)+^(4)+^(7)+j(l) 



127r 



--—+27rp{p~p]^' + 0{^^) 



V2.-K 



(4.19) 



where the last line follows from the compressibility sum rule (3.5) and where we have noted that p as well as p differ 
of respectively and po by terms of order 7^/^. The Legendre transform of usrpm is achieved by combining Eqs. 
( 4.17 ), ( 4.19 ) and by noting that the correlation functions /iq"'' which enter terms of order 7^/^ or higher can safely 
be replaced by the /i*'^"^ since these functions differ by terms of order 7^/^. Our final result for Jrpm is therefore 



27r(p7)^ / d^rh*^^\r) + 



ifn) 



+ 83/2 (^^)7/2 f ^3^-.^/^*(2)(,) _ ^1/2 (^^)7/2 f h*(^Hl,2,3) 



3 



Kexp-K(ri2 +r23 +^34 + r4i) 4 
X ^ + 0(7 ) 

7'12?'23?'34?'41 



?'12?'23 



(P7)^/' lim / dV2d3f3dV4 h<^\n2)h*(^\r3^) x 



(4.20) 



Some comments are in order. 



(i) The expression (4.20) coincides with-|tliat derived by Stell and Lebowitz by means of diagrammatic expansions 
in the frame of the 7-ordering theory.t3~0 

(ii) We have shown that (at least at order 7^/^ and we are enclined to admit this result at any orders) the 
expansion ( [1.20D is independent of the smearing di amet er a. Therefore the choice a ^ could have been made 
from the beginning, which extends the validity of (4.2C) to other reference potentials than hard core repulsions 
(Lennard- Jones potentials for instance). 



V. CONCLUSION 



In this paper we have derived the 7-expansions of the grand-potential and frea energies of the RPM in the frame 
work of the exact field theoretical formalism exposed in the companion paper.EJ In I we had obtained an activity 
expansion of the same quantities. In both cases of either high temperatures (and arbitrary densities) or low fugacities 
(and arbitrary temperatures) the perturbative expansions are obtained by a cumulant expansion of the exact Sine- 
Gordon representation of the grand partition function S/jpm derived in the first paper. Of course one recovers in 
both cases expressions derived years ago by means of .snnKisticated diagrammatic expansions involving more or less 
complicated resummations of infinite series of graphs .E5Ejuj The field theoretical representation of the grand partition 
function of a theoretically important model such as the RPM (of course other models can be considered too) can 
therefore be used to recover most of the results obtained in the past by more traditional methods of classical statistical 
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mechanics. However we have now at our disposal approximation schemes such that the saddle point approximation, 
the loop expansions etc. which have no real counterparts in the conventional theory of liquids. Concerning the RPM, 
a mean field sine-Gordon approxim|ation of S rpm is perhaps the correct starting point for a full understanding of ionic 
criticality wanted by M. E. Fisher .t3 Work in that direction is currently in progress and will be reported elsewhere. 



APPENDIX A: THE HIERARCHY OF URSELL FUNCTIONS 



Let us consider a classical fluid at inverse temperature /3, chemical potential fi in an external potential ip{r). We 
denote by 'E,[ip] the grand partition function and by 2(r) = exp(— /3(^(r) + v)/X^ the activity {v = and A = de 
Broglie thermal wavelepgW. In the absence of the external field (p, the n-body Ursell (or connected) correlation 
functions are defined asllBEJ 



p("'(i,...,n)= n^Wrr 



lnS[(^] 



V3 = 



It can be shownEl that the pc"'' (1, . . . , n) obey the following hierarchy: 



d V 



(Al) 



(A2) 



Defining now the correlation functions . . . , n) = pi" "* (1, we deduce from Eq. (A2) important relations 

which are repeatedly used in the text, i.e. 



p'J d'r, h('HW) = ^ [p'h^'^ (1,2)) - 2^(2) (1,2) 
dVz dVs (1, 2, 3) = {d^p/dv^) - ?,{dp/dv) + 2p . 



(A3a) 
(A3b) 

(A3c) 



APPENDIX B: EXPRESSION OF THE W^jv'S 



The expressions of 11^2 and W3 which are needed in the text are obtained from Eqs. (2.6). The resulting expressions 
are further simplified by the use of the relations (A3) of Appendix After some elementary algebra one finds 



' 2a2 2a 



4! 



pI 



^'^ /dMdVl, 4^^(1,2)0^(1)0^(2), 
L3 



(Bla) 



+ d^n 0^(1) - ^ / d^n m + ^ I d^n 0^(1) 



Po 



4! a 



Pa 



6! 



L / dV-i d'r, A (pg4^)(l,2)) 0^(1)02(2) 



+ 



+ 



48 

4 
48 



d^n d^f, hl^\l,2) 02(1)04(2) 

dVi d^f2 d^r3 h^o\l,2,3) 02(1)02(2)02(3) 



(Bib) 



where Pq and Pq denote respectively the first and second derivatives of the density of the reference hard sphere fluid 
with respect to the chemical potential vq. 
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APPENDIX C: WICK'S THEOREM 



We collect in this Appendix the averages (. . of the moMTOials of </> which are needed in the text. The following 
results have been obtained by making use of Wick's theoremB3'E2lEil 



mfm^)x. 
m?m')x. 
mfm')x, 

mfm^m^)x. 
mYm?m^m^)x. 



= ^r(O) 

- 3A2(0) 

= 15X^(0) 

= 105X^(0) 

= 2^2(12) +^2(0) 

^ 12X2(12)A^(0) + 3X^(0) 

= 24X^(12) + 72X2(12)A2(0) + 9X^(0) 

^ X3(0) + 8X^(12)X^(23)X^(13) 

+ 2X,(0)[X2(12) + ^2(23) + ^2(13)] 

= 3A^(0) + 48X^(12)X^(23)X^(13)X^(0) + 24X2(12)X2(13) 
+ 12X2(0)[X2(12) + X2(13)] + 6X2(23)X2(0) 
= ^.'(0) 

+ 2X2(0)[X2(12) + X2(13) + X2(14) + X2(23) + ^2(24) + ^2(34)] 

+ 8X^(0)[X^(12)A:^(23)A:^(31) + Xr{l2)Xr{24)Xr{'il) 

+ X,(13)X,(34)X,(41) + X,(23)X,(34)X,(42)] 

+ 4[X2(12)X2(34) + X2(13)X2(24) + X^^{U)X^^{2'i)] 

+ 16[X^(12)X^(23)X^(34)X^(41) + A^(12)A:^(24)A:^(43)X^(31) + 

+ X^(13)X^(32)X^(24)X^(41)] 



(CI) 
(C2) 
(C3) 
(C4) 
(C5) 
(C6) 
(C7) 

(C8) 

(C9) 



(CIO) 



APPENDIX D: COMPUTATION OF {Wi)x^ 



The cumulant (Wf )xt ^ easily obtained from ( Bla ) and reads 



X' 



X 



r,'2 



4a2 
(4!)2 



= / '^'^^'^''^^ {{4>\l)<t>^{2))x. - {<t>^{l))xA<t>^{2))x^ 



d'nd'r2 ((<^^(l)0^(2))x. - (04(l))x.(04(2));,J 



X3 = g / d\-,d\-,d'r,d'u 4^^(1,2)4^^(3,4) ((,/>2(l)02(2)02(3)^2(4)^^^ 



.3 _ P|y dV-idVidV-gdV-i 4")' 

b^l)q^^{2))xA<P'{3)^'mx.) , 



4!a 



dV-ldVl ((02(1)^4(2))^^ _ (02(1))^^ (04(2)) 



8a 



- ^PoPo / ^3^ .3-. -3^„ . (2) 



d'nd^r^d'fs C(2,3) ((02(1)^2(2)02(3))^^ _ (02(i))^^(02(2)02(3))^j 



6 _ PO 



d'nd^r^d'r, 4')(2,3) ((04(1)02(2)02(3))^^ „ (04(1))^^ (0^(2)02(3));,^ 



4x4! 

Each of the X* are now computed by making use of Wick's theorem. It follows from (CI) and (C5) that 

L-^X' = -^/(2) , 
2a^ 



(Dla) 

(Dlb) 
(Die) 



(Did) 
(Die) 
(Dlf) 

(D2) 
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where 



= J d^rX^r) . 
Similarly, making use of Eq. (C2) and (C7) one finds for 

8 24 



with 



= J cfrXtir) 



(D3) 



(D4) 



(D5) 



The term X^ is considerably more complicated. It follows from Eqs. (C5) and (CIO) that it can be splitted into four 
pieces 



(D6) 



We have defined X'^^°''> as 



(D7) 



where we ma de u se of the compressibility relation ( |3.5| ) to obtain the last line of Eq. (D7). The second contribution 
X^C*) in Eq. (Q reads as 



= -plXr{Q) / d■'r^dPr2d'T:id'r^ /ij,^^ (1, 2)/i^^^(3, 4)X,(1, 2)X,(2, 3)X,(3, 1) 



^r(O) (p; - Po) 1 d'V-'4^)(r)X.(r)X;(2)(,) , 
where we have introduced the auto-convolution 

= j d'^T^ Xr{r)Xr{r-r) . 

The third contribution X^'^'^^ is given by 

= dVldV-^d^fgdVl 4^)(1,2)4^)(3,4)X,2(1,3)X,2(2,4) 

which can be recast under the simple form 

x3w=^i3 j d^f(^h^^Kx^y , 

where * X^ denotes the convolution of the functions h^\r) and X^{r). Finally the last contribution X'^f'') reads 
as 

Xm ^dj d^nd'f2d''r3d'ri 4'^(1,2)4'^(3,4) 2)X,(2, 3)X,(3, 4)X,(4, 1) 

+ X,(1,2)X,(2,4)X,(4,3)X,(3,1)+X,(1,3)X,(3,2)X,(2,4)X,(4,1)] , 
which can be rewritten as 



(D8) 



(D9) 



(DIO) 



(Dll) 



(D12) 



11 



3(d) ^P0^3 I ^3^{{h(^) 



(r) 



(ffid^r2(fr3d^ri (1, 2)/i^^^ (3, 4)X^(1, 3)X^(3, 2)X^(2, 4)X^(4, 1) , 



,(2) 



(D13) 



(2) 

where, once again, we have represented the convolution of the two functions Hq {r)XT{r) and Xr{r) by the symbol 



Let us quote now the expressions of the remaining without further comments 
L-^X" 



L-^X"" 



1 , 

2a 



Pi 



^X.(O) j d'fhi^\r)XAr)X<'Hr) 
Collecting the expressions of the X* one finds finally that 
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i=l 



w'2 = A2(a) / d'rh^^\r)Xr{r)X<'Hr) 



w,^Pf ld^rh^^\r){Xl . Xl){r) 



pI 



^4 = ^ * X,-) (r) 



^^^y J d'r {(h^^^X,) * Xr) (r) 

= ^ J d^nd^r^d^d^'n 4''(1,2)4'^(3,4)X,(1,3)X,(3,2)X,(2,4)X,(4,1) 



W7 



24 



^i(«) = Ipo'^l + Oia) 



1 



^2(a) = -|p§Po'«o + 0{a) 



(D14a) 
(D14b) 

(D14c) 

(D15a) 

(D15b) 
(D15c) 

(D15d) 

(D15e) 

(D15f) 

(D15g) 

(D15h) 



with complicated expressions for Ai{a) and ^2(0). However both functions are regular in the limit a — > 0; more 
precisely one has 



(D16a) 
(D16b) 



We now proceeds in two steps. We first examine the limit a ^ of each of the term in the r.h.s of Eq. (D15a) and 
then study the limit 7 ^ of the resulting expression to decide whether or not it contributes to "j'^ {W2) at 
order C'(7''/2). In other words, we seek the terms of limQ^o(W^2 which are singular in the limit 7 — > and discard 
the regular ones. 

Since in the li mit a the potential Xr{r) tends to the Yukawa potential y{r) = exp(— Kor)/r one has, taking 
into account Eq. (D16a) 



P0..2 /■ j3-exp(-2Kor) _ 7Tp„ 



lim wi = / d'^r 



(D17) 



and therefore the term w\ can be discarded. Hag aEi has shown that y*^'^\r) = 27rexp(— Kor)/Ko from which it follows 
that 



12 



lim W2 

a-*0 



r 



(D18) 



Since h^\r) is a short range function of r the integral in Eq. ( D18| ) converges and the term W2 can also be discarded. 
In the limit a — s- the third term 1^3 reads 



1- Po 

iim W3 = — 

a^O 4 



y*(2) ^ y*{2) 



(r). 



(D19) 



It turns out that the function K{r) = y*^^^ * y*^'^\r) has also been studied by HagaEj who proved that K{r) = L(r)/r 
where L(r) is an entire function of Kpr, from which we conclude that lima^o is regular in the limit 7^0 and can 
thus be safely discarded. 

In order to study limQ^g W4 we rewrite the integral in Fourier space and obtain 



Pt 

lim WA = — 

a^O 8 



(2^ 



U;^\q) ~[TF{y^){q)\ 



The Fourier transform of y'^{r) has a simple analytical expression, also given by Hag ai, h 

TF{y^){q) 

and therefore 



47r q 

— arctan 

q 2ko 



(D20) 



(D21) 



lim lim W4 

7^0 a^O 



4 4 

Pa 



d^'q 



(27r)3 



(D22) 



Since the integral (D22) is convergent the term W4 can, as the previous ones, be discarded. 
In order to study the limit a — s- of the term W5 we have to consider the integral 



% {r)y{r) * y{r) 



(D23) 



If we interpret z(r) = {r)y{r) as a spherically symmetric distribution of Yukawa charges, then v{r) = z * y{r) is 
t he p otential created by this distribution. Obviously v{r) ~exp(— Kor)/r for r 00. Therefore the integrand of eq. 
( pl)23| ) behaves as exp(— 2Kor)/r^ for r 00. The integral ( D23 ) diverges thus as I/kq in the limit 7 — > 0. It follows 
from this remark that 



lim 7 W5 

a— *0 



jpoir^' lini 



d'nd^r^d'r, h^^\l, 2)h'P (3, 4) "° ^"P^-""^'^^^ + + + ^^1)) 

7'127-237'34?'41 



We come now to the term wg . Taking advantage of the translational invariance we have 



p4 

lim lim wq — — 

Ko^Oa^O 4 



,(2)/ 



1111 

f 13 f 32 ^'24 rii 

(2), 



(D24) 



(D25) 



which is a convergent integral due to the short range of the pair correlation (r) . Therefore wq does not contribute 
to ujRPM at order 7''/^. 

The last term to examine is w-j which is a tricky one, since the integral /^^^ div erge s at short distances in the limit 
a — > 0. In order to study carefully this divergence we insert in the expression (|3.9[) of XT-ir) which yields 



hm = ^ - + 16k^ (7b - 1 + In (aKo)) + 0{a) 



(D26) 



where is Euler's constant. Note that the divergence of ^'^wi occurs at order 7^ and that we do not need to take 
it into account at the order of our calculation (i.e. 7^^^). However we have checked that this divergence for a ^ 
originates from the fact that, till now, we have ignored contributions coming from the cumulant ^^{W/^) x^- Taking 
carefully account of one of the terms of this cumulant leads to consider the integral 
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/(4) ^ 1 ,3,. exp(-4.or) ^(.) ^ ^^^^^ ^^^27) 

instead of I^. Recall that 5^^^(r) = h^o\r) + 1 is zero in the core which makes / ^ ' convergent in the limit a ^ 
and lim^^o linia-»o I perfectly well defined. However the corresponding term j'^lo-j contributes to lvrpm at order 
7* and can thus be discarded. 

The conclusion of this lengthy discussion is that 

Ko exp{-Ko{ri2 + r23 + 7-34 + r4i)) „ , 4 



r-i2r-23r34r4i 



+ 0(7^) . (D28) 
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